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Abstract 

The Navier-Stokes equation for incompressible liquid is considered in the limit of 
infinitely large Reynolds number. It is assumed that the flow instability leads to gener- 
ation of steady-state large-scale pulsations. The excitation and evolution of the small- 
scale turbulence is investigated. It is shown that the developed small-scale pulsations 
are intermittent. The maximal amplitude of the vorticity fluctuations is reached along 
the vortex filaments. Basing on the obtained solution, the pair correlation function in 
the limit r — > is calculated. It is shown that the function obeys the Kolmogorov law 

r 2/3^ 

1 Introduction 

In the turbulent flow, in addition to the average velocity of the flow, random velocity pulsa- 
tions are excited. These pulsations could be presented as a sum of different scales random 
movements. The large-scale pulsations with scale L of the same order as the characteristic 
parameters of the flow play a leading role (For example, in the tube of radius R the scale is 
L ~ R/5). Large-scale pulsations have the highest amplitudes. 

The small-scale pulsations with scales I « L are excited also. They have much smaller 
velocity amplitudes, and they could be considered as a fine structure set to the main large- 
scale movement. The small-scale pulsations contain only a small part of the whole turbulent 
kinetic energy (see Landau, Lifchitz [T], Monin, Yaglom [2]). 

If the viscosity of the liquid v is small enough, and the Reynolds number R u correspond- 
ingly is large, then the spectrum of the small-scale pulsations becomes very wide. This type 
of turbulence is called developed. Let Ao be the maximal scale where the viscosity is still sig- 
nificant; then the range of scales Ao << / << L is called the inertial interval. The pulsations 
developed inside the inertial interval in different scales are determined by nonlinear processes 
only, since the viscosity v is negligible. Therefore, it is possible to study the inertial interval 
of the turbulence in the limit v — > , A ^ , R u — ► oo (Frisch [3]). 

In the turbulent flow velocities are random. So, the correlation functions could be used 
to describe them. Let us consider the isotropic turbulence. The pair correlation function 

ir(r) = <[v(p)-v(p + r)] 2 >, 
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determines the relation between the values of velocity in two near points p and p + r. Since 
the turbulent pulsations are isotropic, the correlation function depends on the distance r 
between the points only. The pair correlation function measured in numerous experiments 
has the universal form: 



The distance is restricted by the condition r << L, i.e. the experimental result (CQ) refers 
to the inertial interval only. The Fourier-transform S(k) of the correlation function (CQ) has 
been also investigated experimentally. These investigations give 



where k is the wave vector. The spectrum (j2J) is called the five-thirds law. The limit r — > 
corresponds to k — > oo. In the developed turbulence the five-thirds law is observed inside a 
wide range of wave numbers - up to three-four orders of magnitude [2J, [3]. 

We emphasize that the experimental measurements of the correlation function in the 
small-scale region were initiated by the theoretical predictions. A.N. Kolmogorov in his 
fundamental works in 1941 [I] derived the expressions O, (J2J) for the velocity correlation 
function in the homogeneous and isotropic turbulence. □ The Kolmogorov's theory is phe- 
nomenological. Its basic conception is the uniform dissipation of energy in the turbulent 
liquid. There is a stationary flux of energy in the Fourier space: the energy is generated 
in large-scale pulsations, and flows uniformly through the whole inertial interval of scales. 
In this process, the flux of the energy is conserved. The dissipation occurs only outside the 
inertial interval, at the smallest " dissipative" scales < Aq- Relaying on this physical model, 
using the relations of similarity and dimensions and the general properties of hydrodynamics 
equations, the correlation function was found. 

This fundamental Kolmogorov's result was later confirmed in numerous experiments. It 
stimulated a huge amount of theoretical, mathematical and (in the recent time) numeri- 
cal investigations. In these works the theory of turbulence was widely developed (see the 
monographs [2], [3], [6] - [9] and literature therein). Recently, the methods of field theory and 
solid-state physics have been used [TD],[TT]. However, the attempts to obtain the expression 
for the correlation function directly from the Navier-Stokes equation without any additional 
assumptions have not been successful up to now (see [3] for more details). 

Another approach to the problem is based on the physical ideas of the leading role of 
singularity in the small-scale structures of developed turbulence [3], [13]. However, despite 
significant efforts in this direction, neither the correlation function has been derived from 
the Navier-Stokes equation, nor even the existence of singular solutions has been proved. 

Thus, the problem of derivation of the fundamental Kolmogorov's result directly from 
hydrodynamic equations is not solved yet (see the monographs [2J, [3]). 

In this paper we propose a new approach to the problem. It allows to find the structure 
of the small-scale turbulence and the pair correlation function. 

From the hydrodynamic equations written in the Lagrangian reference frame, we derive 
the equation describing the joint probability density of vorticity u) = V x v and its time 

1 The law ([5)1. which is the direct consequence of (TT]), was written in an explicit form in the papers by 
A.M. Obukhov [5]. 





S(k) = C f k' 5/3 , 



(2) 
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derivative. We show that moments of the vorticity distribution grow unrestrictedly in time. 
Then we find an asymptotic solution at infinitely large time. Basing on it, we obtain the 
spatial distribution of the vorticity where it is large. These are vortex filaments. They give 
the main contribution to the pair correlation function. 
The paper is organized as follows. 

In Section 2 the equations of motion of incompressible liquid are considered. Their 
decomposition in the vicinity of trajectory of an arbitrary lagrangian particle is written. It 
is shown that local vorticity growth is determined by anisotropic part of large-scale pulsations 
of pressure. 

In Section 3, supposing the randomness of large-scale pulsations of pressure, the equation 
for probability density of vorticity and its time derivative is obtained. We show that even 
moments of the vorticity distribution grow exponentially, the higher moments growing faster 
than the lower ones. This is the manifestation of intermittency of hydrodynamic turbulence 
in small scales. 

In Section 4 we find the large time asymptotic solution for the joint probability density 
of vorticity and its time derivative. 

In Section 5, on the ground of the obtained asymptotic solution, the spatial structures 
contributing mainly to the asymptote of the probability density are investigated. We show 
that these are the vortex filament structures which determine the pair correlation function 
of turbulent pulsations p]) as r — > 0. 

In Conclusion we formulate and discuss the main results of the paper. 



2 The statement of the problem 

Let us consider the Navier-Stokes equation for incompressible liquid. It is known that at the 
scales larger than the viscous scale Ao it takes the form of the Euler equation: 

^ + (v-V)v + ^ = 0; V-v = (3) 

at p 

Here v is the velocity of the flow, p is the pressure. The density p is taken unity below. The 
second equation expresses the incompressibility of the liquid. The equations ([3]) describe 
the processes on the scales inside the inertial interval (see [3]). From ([3]) one can find the 
relation connecting the pressure with the flow velocity: 

— Ap = V{Vj ■ VjVi 

To investigate the local properties of the turbulent flow we pass on to a coordinate system 
co-moving to some element of the liquid with coordinates $(t): 



r' = r - £(t) , v = v - i; £ = -Vp 



Dot means the time derivative. After such change of variables the equations d3J) take the 
form 

dv' 

— + (v- V)v' + VP = 0, VP = Vp + £, W = (4) 
at 
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Since the reference frame is chosen to be co-moving, at the point r = we have 



VP(r = 0) 







v'fr' 



0) = 



Expanding the velocity v' and the pressure P into a Taylor series in the vicinity of the 
co-moving point we have the main term: 
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P = ^pikr H r' k , p ik 



VNkP 



(5) 



(6) 



Here the tensor 



dr'i 



is decomposed into a sum of symmetric and antisymmetric 



\tikjXk parts. Note that ba = since Vv' = 0; pij is symmetric. It is easy to check that the 
vector Xi defined by the asymmetric part of dv'Jdr'i is the vorticity (uj = V x v) of the flux 
at the point 







Xi 



Combining (jSl),® and (J4]), we obtain 



1 2 \ 1 

bij ~ {Xi'J'j X ~\~ bikbkj -\- Pij J -\- — (Sikjik £ jkn%kbin £-ikn%kbjn) 







The first term in the square brackets is symmetric, the second term is antisymmetric. Since 
the equation should hold for any r', both terms are equal to zero. Actually, multiplying the 
expression in the brackets by ey„ and taking ba = into account, we find for Xi and b^: 



bij + ^(xiXj - x 2 5ij) + b ik b kj + = 



(7) 



b n k%k (8) 

Taking the time derivative of (jSD, we obtain finally the equations set for the three components 



Xi . 



i>n Pnk-Ek (9) 

Let us clarify now the physical meaning of the symmetric part of velocity b^. Namely, 
let us express it in terms of space distribution of the vorticity u>(r). We shall see that in 
the completely isotropic flow b^ = 0. Hence in accordance with (jSJ), the vorticity at the 
lagrangian point does not change. In the real flow there are two regions where the isotropy 
may be broken: either local, on account of small-scale pulsations of pressure in the vicinity 
of the point under consideration; or global, the remote areas close to the boundary of the 
system, at the scales of the order R. We shall show that is determined just by the global 
break of isotropy. 

Since V • v = 0, there exists a vector potential A: 



V x A, V - A = 0. 
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Then 

AA = -u>. (10) 

To separate the singularity accurately, let us expand A(r) and u>(r) into a series on spherical 
harmonics: 

oo m=l oo m=l 

A = E A im (r)y /m (0, v?) , u> = w lTn (r)^ m (0,¥>). 

Z=l m=—l 1=1 m=—l 

The solution of the Poisson equation (|TU|) is: 

A im (r) = ^ u^nY^dn + ^-j- y ^^(rOrJ-'dri (11) 

The integration limits are chosen to provide the convergence of the integrals at r —>■ and 
r — > oo. Note that for analytic function u>z m (r) oc as r — > 0. 

To evaluate by, we need to determine the limit VfcU, = SynVfeVjAn as r — > 0. Only 
quadratic (in the coordinates r*) part of A contributes in it. This quadratic part consists of 
two terms proportional to the zeroth A o and the second A 2m spherical harmonics. Hence, 
we are interested in the two harmonics only. 

The zeroth harmonic A o gives the local contribution corresponding to the antisymmetric 
part of the velocity tensor: x = u;oo| r= o = ~~ A A o | r= o- However, in the symmetric tensor 
b^ its quadratic component A o oc r 2 is cancelled: it is just the fact that gives b^ = in 
isotropic medium. So, only the second harmonic A2 m remains. Since u>2™(r) ~ f 2 as r — > 0, 
we see that the first integral in ffTT]) behaves like r 4 , and the second one - like r 2 as r — > 0. 
Hence, the contribution of the "local" item is negligibly small, and the symmetric part of 
the velocity tensor by is determined by "global", large-scale properties of the whole flow. 
Returning to the rectangular coordinates and taking the derivative, we obtain 

K = e jnk J (** - 3^) dr> + (i ~ j) 

According to our analysis, the integrand has no singularity at r = 0; the integral accumulates 
at the scales of the order of R, where the isotropy breaks. 

The analogous argumentation shows that the pulsations of pressure (see (JTj)) could 
also be presented as a sum of local and large-scale pulsations; the local part of the tensor 
is XjXfc — SikX 2 . From it follows that this tensor does not affect the local vorticity of the 
flow. Hence, the local dynamics of vorticity depends on the large-scale pulsations of the 
pressure p n k only. 

Thus we obtain the first main property of the turbulent flow: the local vorticity along the 
streamline in homogeneous and isotropic flow is determined by anisotropic part of large-scale 
pulsations of the pressure. 



2 Note that in two-dimensional flow such a division into local and large-scale components is not possible. 
The zeroth cylindrical harmonic A diverges logarithmically, and as a result the "local" component should 
influence on the large-scale component. 
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3 Probability density equation 



Since we interested in statistical properties of the flow, let us introduce the probabilistic 
description. We consider now the vorticity u>(£) as a random quantity. Its change still obeys 
to (El). Instead of one equation of the second order, let us consider a system of two first-order 
equations: 

Xi = yi yi = -pijXj (12) 
Here = Ui and i/i = u>i. We introduce a joint probability density 

f(t, x, y) =< 5(x - x(t))5(y - y(t))) > . (13) 

Here x(t),y(i) are the solutions of ([12"]) at the given realization of pij and initial conditions; 
the average is taken over the ensemble of all possible realizations. 

The aim of the paper is to study a steady-state turbulent flow in the inertial interval of 
scales, i.e. at scales I and time t satisfying to the conditions 

I « L, t » r c (14) 

Here L and r c are the characteristic space and time correlation scales of the large-scale 
vortices. These large-scale vortices depend on the specific geometry of the installation and 
on the boundary conditions. According to the experimental data [Hj, the large-scale velocity 
pulsations are random and Gaussian. Thus, in the equations (Q or ffl2l) for the local vorticity, 
the matrix Pijit) describing the large-scale fluctuations of pressure could be taken Gaussian 
and, because of (1T4"1) . delta-correlated in time. These propositions would be discussed below. 

Note that, as it follows from ([9]), the "random" behavior of vorticity (or velocity) is 
caused by the randomness of the large-scale flow and the corresponding matrix pi k {t). 

The Gaussian random process is described by a pair correlation function 



< Pij(t)Pki(t') >= D ijkl S(t - f) (15) 
Using ( fl2i) and taking time derivative of the probability density function, we obtain 

df df d 

~dt + Vk dx^ = Xp dy~ k ( pkp6 (* ~ x W) 5 (y ~ y(*))> ( 16 ) 

Let i?(x, y,p) be a functional of p. To find the correlation function (p kp R(x., y, p)) we use 
the standard averaging technics for delta-correlated random process (see the monograph by 
Klyackin, [15]): 

(z k R[z\) = W dt' {z k (t)z k ,(t')) /6R[z ' t] 
hi ^ 



k 

Taking ffl5|) into account, we get: 



5zy{t') 



<p fcpJ R(x,y,p)> = ]TZW ( S f [ * y ;J ] ) ( 17 ) 
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To evaluate the variational derivative (TITl) . we use the equations of motion ( fl2l) : it follows 



Sp k 'p'(t') 



-Skk'Xp'it) 



Sx k (t) 



t=f 



5p k y(t') 



t=t' 



Combining this with ffTBT) . we obtain the Fokker-Planck equation for the function /(t,x, y): 



dt ^ k dx k 



DijklXjXl 



d 2 f 

dyidy k 



The matrix Pij(t) is symmetric ([6]). Hence, in the homogeneous and isotropic medium the 
general form of the matrix Dij k i is 



D ijk i = D5(t- t') (8 ik Sji + SiiS jk + TSijSki) 



(19) 



The constants D and T depend on the large-scale flow. In addition to the isotropy and 
homogeneity, it is natural to suppose statistical independence of different components pa-. 
In this case one has T = 0. However, the values D and F appear to be unimportant. The 
parameter D in the equation vanishes as a result of time normalization. As we shall see 
below, the resulting properties of the turbulence depend only weakly on the parameter T 
(the positive defmiteness leads to a restriction T > —2). 
Substituting ( fl9l into ( |T8l) . we obtain finally: 



df , df 
at oxk 



2 d 2 f 



x 



<9y 2 



+ 7 \Xk 



d_ 

: dy k 



f 



(20) 



Here 7 = 1 + T, time t is normalized by D 1 ^ 3 . The value D^ 1 ^ 3 is the characteristic time of 
probability density change. As it was shown above, in the completely isotropic turbulence 
D — 0. Taking into account small anisotropy, we have D^ 1 / 3 » r c . This allows to use the 
delta-correlation approximation in derivation of the equation (|20|) . 
We now itemize the main properties of the equation (|20|) : 

1. All momenta of values Xk and yj of the order n are connected by a system of the first 
order linear differential equations. 

2. Even momenta grow exponentially. Independently of the initial conditions, the func- 
tion / at large values t depends on the modules x, y and the cosine of the angle between the 
vectors p = (x, y)/xy only. 

3. The higher even momenta grow faster than the lower ones. 

To illustrate these statements, consider the momenta of the second and the fourth order. 
Integrating the equation (|2"U1) in x and in y, we obtain for the second-order momenta: 



d 



dl 



< XiVj + xjyi >= 2 < y i y j > 



(21) 



dt 



< ViVj > 



2§ij < x 2 > +27 < XiXj > 
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Let us consider the invariant momenta of the second order, i.e. < x 2 > , < y 2 > and 
< x • y >. Their evolution is determined by characteristic equation 



A 3 -4r-16 = 0. 

Asymptotically as t — > oo we obtain 

< x 2 >oc< y 2 >oc< x ■ y >oc exp(A 2 t) 

Here A 2 = (16 + 4F) 1 / 3 . 

For the other momenta (i ^ j) one has 

a 2 = 4 + 4r 

We see that A 2 > A 2 . Hence, at large time the invariant momenta are much larger than 
the others. In other words, the probability density function at large time depends on three 
variables x,y, fi = Xiyij [xy] only. 

The characteristic equation for invariant momenta of the fourth order takes the form: 

A 6 - (84r + 244)A 3 - 1280 = 0. 

For example, for < x 4 > one has < x 4 >oc exp(A 4 t), where 



A, 



1/2 lV3 



((84r + 224) 2 + 5120) ' + (84T + 224) 



One can check up that for the values T > —0.9 holds A 4 > 2A 2 . Hence, as t — > 00 one has 

< x 4 >>< x 2 > 2 . 

The obtained relations demonstrate that the higher momenta of the vorticity module 
grow exponentially in time. This property is called intermittency. It reveals itself in the 
instability of small-scale flow. Physically this instability means that under the influence of 
large-scale random pulsations a drop of incompressible liquid stretches out. This leads to 
generation of vortex filaments. These filaments provide the basis of the small-scale turbulence 
(a simple physical example demonstrating the process of filament growing is considered in 
Appendix). 

It will be shown later that the domain of parameters y 3> x plays an especially important 
role. In this domain the suggestion of gaussian random process is not needed. Actually, from 
f[T2l it follows that the change of y during the correlation time is Ay ~ x. Hence, 

Ay x . . 

— y - « - < 1 if y^>x 22 

y y 

In this case the fluctuations of the probability function are very small: 5f / >, and one 
can obtain the equation (1201) using the perturbation theory. As a result, the equation has 
the Fokker-Planck form. 
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4 Asymptotic form of the probability density function 



As it was shown in the previous section, the probability density function (1201) at large time 
depends on three variables only: /(x, y) = f(x,y,fi), where \i = (x.,y)/xy. Besides, the 
equation ( l20l) and the initial conditions to the probability function allows integrating over 
three other variables. 

As a result, the equation (1201) takes the form: 



£ + ™ ^ - S£ (4) 4U {i -<) (23) 

The function / must satisfy the normalization condition J fdxdy = J fx 2 y 2 dxdydfi = 1, 
and two conditions of zero flux from the boundaries x = and y = 0. Let us specify the 
meaning of these conditions. For that we return to the equation (J2DJ). It has the divergent 
form df /dt = V Q J Q , a = 1..6. The flux density J a in a 6-dimensional space (x, y) is 

2 df , . (_ df 



-yf,x — +7X x 

oy V <?y 

The no-flow boundary condition at y = means that the integral of J over the 5-dimensional 
surface |y| = e vanishes as e — > 0. After integrating in all angles dfl x dfl y = 4n ■ 2ndfi, in 
terms of variables x, y, fi we have 



2 9f , . 
x — + 7X ( x 



^r- I }--x 2 dxdtt x y 2 dtt v = 8rt 2 f (l + ^a 2 )^x 2 y 2 dxdu — ► 0(24) 

dy/J y J d y y _> 



Similarly, the no-flow condition at x = leads to 

fuyfx 2 y 2 dfidy — > (25) 
x -> 

The expressions (12~4"|) and (1231) are the boundary conditions for the equation (1251) . 

We now return to (T23|) . We search for a stationary solution as t — > oo. Choosing a new 
variables x, z = y 3 /(3x 3 ), fi let us present the function f(x,y,fi) in the form 

/(x, J/, A*) = x ~ 2x ~ aF ( z i M a) (26) 

Q: 

The set of eigenvalues a is to be found by solution of (12"51) with boundary conditions 
The equation (123)) then takes the form 



d 2 F (A \ dF a _ 1 9 r/1 2N _, 19 
* fl? + U + "*J & + 3 ^ " 3 ^ " " ^ + 9^ 



'1 



(27) 
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+7 



a | i-/i 2 d y p | 2 / g 1 - a* 2 a 

cte 32; <9/x/ 3 \ dz 3z d\i 



Integrate (|27|) on variable /x and define functions (j,(z) (J> 2 {z): 
Jt(z) 



J^fiFdfi Ft — , J^fi 2 Fdfi 2F 2 1 
■ + - 



/x 2 (zg 



f28l 



Here -Ffe(^) are coefficients in Legendre expansion of the function F. The equation ([2"T]) takes 
the form 



2(1 + jh 2 )F 0zz + [ \ +JIz + 7 ^/i 2 2 + T/i2 1 | | F 0z 

-o 7-5 3]? - 1 
+ I -/x + z/x^ + 7 I z/x 2 22 + -/i 2 2 + 



3' 



92 



Fn = 



Substituting F (z) = w(z) exp ( — f* \ , we get 



(29) 



, . 4 7/i 2 — 1 — \ /a 

(1 + ^)zw zz +[- + + 2 7 ,/x 2 , -nz]w z + + ^_ 



/xw (30) 



'l-"" + 3 ' i * + 9z l + 7 ^ 3(l+ 7 ^ 



=^ w = 



n 00 

The solutions of fl30|) could be presented as a series □ to = z s c n zH - It converges on the 



n=0 



domain < z < oo (if 7u(z) and /x 2 (2:) have no singularity). 

In order to find s and c n , let us consider the asymptote 2 — > 00 @. Expanding (1291) into 
Legendre series, we get 



4 1 r 

+ --TOz + g 



ft 



zF l2 + -F x 







+ + (4 + 1 ) ( F ° + l F >) + l F " ~ TE F > = S Fl 



Z ^m + q -^m + ( Z o q I 1 

3 \ 3 / V 2m- 1 



m+1 

'm— 1 "I - ^ m+1 

2m+3 



1 / m(m+l) m(m+l) \ m(m+l) 

"I - " I m— 1 F m+1 I « 



3 V 2m-l 



2m+3 



9z 



7 / m 2 (m-l)(m-2) 2m(m+l)-l (m+1) 2 (m+2)(m+3) 
+ — I -Tm-2 — m{m + L) t m H -Tm+2 



92; V (2to-1)(2to-3) 



(2m-l)(2m+3) 



(2m+3)(2m+5) 



3 Zero is a regular critical point of (fHUj) . since 7Z(z) and ^ 2 (z) are unbounded, and 1 + 7 > 0. 

4 As it has been mentioned in the end of previous section, in this limit the equation for probability function 
has the Fokker-Planck form independently of statistical properties of the large-scale random process 
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, m(m-l) „ 2m(m+l)-l „ (m+2)(m+l) „ 

+7^ I ^ m—1 + ' m + £ m-\-1 

1 (2m-l)(2m-3) (2m-l)(2m+3) (2m+3)(2m+5) 



7 / m(m— l)(2m-5) , 8m(m+l)-4 , (m+l)(m+2)(2m+7) 

+ 77 I m—2 ' m ' ■ 

3 \ (2m-l)(2m-3) (2m-l)(2m+3) (2m+3)(2m+5) 

As z — > oo, we neglect the terms proportional to (-^f). The resulting equations set has the 
solution 

F m = {2m+l)F {z) (31) 

where Fo(^) satisfies (1291) for /J = /i 2 = 1. (Actually, the coefficients (13T1) are the Legen- 
dre coefficients of the function 26(1 — fi)F (z; a).) Combining (I3TT) with the definition of 
Ti{z)^(z) (EHD, we get 

?*(*) = = i-oQ 

The equation fl30|) takes the form 

(1 + i)zw zz + \- + 2-f-z)w z + ( — - ) u; + 27 | — - — — — - ) w = 

1 J V3 7 V3 3(1 + 7) J J \9z (I+7); 

This is equivalent to Kummer degenerate hypergeometric equation [IB] . The solutions of 
this equation are 

= ^-It^M fa, 6; ( 32 ) 



w 2 (z) = Z - 1/3 M (l + a - b, 2 - b- 



where M is the Kummer function 



, r , , >n -, «z a(a + 1) -z 2 



the parameters a and 6 are 



4-a , 22+7 
= -- 



3 3 1 + 7 

We have found the general solution. Let us now check the boundary conditions 
). The solution w 2 gives F ~ z~ 1//3 as 2 — > 0. The correspondent J y 2 ^dfi does not 
vanish as y — > 0. This contradicts to the boundary condition (1241) . Hence, the solution of 
our problem is u>i, since it satisfies ( |24"|) . 

The Kummer functions behave like M(a, b, z) ~ e 2 z a_b as z — > 00 if a is not negative 
integer. Therefore, for the corresponding values a = 4 — 3a and for small x we have F(z) ~ 
z _a//3 , /x 2 ~ y~ a . This means that the no-flux condition on the boundary x = ( 1251) is 
not satisfied. Hence, to satisfy both boundary conditions ( 1241) and (J25j) one should take the 
solution w\ with the values a that correspond to "discrete" spectrum □ 

a = 4 + 3n, n = 1,2,3.... (33) 



= is excluded since it does not satisfy the normalization condition for /. 
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For these values a the series M contains a finite number of terms, the leading term being 
~ z n . 

The solution 0321) together with (1331) and (1311) gives asymptotic behavior of the function 
/ as t — ► oo. We stress that, according to (1261) . fl33l) . the full probability density function 
f(x,y,fj) in the leading asymptotic term behaves like 

f(x,y,n) ~ x~ 9 F(z,[i; 7) 

as x — > oo. Below, we will need integrability of the function F only. 



5 Spatial distribution of vorticity. Singularity of vor- 
ticity and pair correlation function. 

In the previous section we found an asymptotic solution for probability distribution of u>, u>. 
It is important that the solution has power fall as x — > oo. This means a significant prob- 
ability of large-amplitude fluctuations of \u>\. This is the manifestation of intermittency in 
the turbulence: in some spatial domains the value |u;| is much larger than its average. The 
question is how the flow in this domains should look like to provide the obtained asymptote. 

For this purpose, let us define the probability density function for the module of vorticity 
based on the combined probability density of u>, ci> (fT3"j) . (}2"U]) : 



P{x,t) = J f(t,x,y)x 2 dydv=(6{x-\x(t)\)) (34) 

On the other hand, we independently define a probability density Pi(x, t) as space average 
of some realization of the turbulent flow: 

Pi(x, t) = ^ j 5(x- X{t, r))dr. (35) 

Here X(t, r) = \u)(t, r)| is the vorticity module at time t and at the point r, V is the volume 
of space occupied by the flow. 

The first expression for the probability density P(x, t) is the ensemble average along a 
trajectory of liquid particle and the second expression for Pi(x, t) is the space average. 

Owing to ergodicity (i.e. the equality of ensemble and space averages) we obtain 

P(x,t) = Pi(x,t) (36) 

If the function P(x, t) is known, then it is possible to derive space distribution of the 
vorticity module X(t, r) using fl36l) and fl35|) . Since we are interested in possible singularities 
and their surrounding, let us consider the limit t — > oo and next x — ■> oo. 

Suppose that the singularity is reached at some surface. Taking the point of origin on 
the surface and the axis z perpendicular to it, we find 



If 1 

Pi(x, t) = — / 5 (x — X(t, z)) dodz = — — 
V I \X' 



X(t,z)=x 
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Here da is the element of the surface area. The simplest example of surface where the 
vorticity grows unrestrictedly large is the tangential break of flow velocity. Indeed, in the 
case the velocity is V Q on one side of the contact surface and — V on the other side. Hence, 
the vorticity is concentrated on the surface. 

Let us now consider the most interesting case: the maximum of X(t, r) is reached along 
a vortex line. Then, choosing the cylindrical variables z, r, <p with z axis oriented along the 
line, we obtain 



Pi(x,t) — — J 5 (x — X(t,r))rdrd(pdz 



\XL 



X(t,r)=x 

In the case of point-like maximum, using spherical coordinates r, 6, 0, we would obtain 

,2 



Pi(x,t) = — J 5 (x — X(t, r)) r 2 dr sin 



XL 



X(t,r)=x 



Taking the limit t — ► oo, with account of (1361) we find the spatial distribution of vorticity 
module in the vicinity of singularity: 

X'(z)P(X) = 1 — singular surface 

X'(r±)P(X) = r±_ — singular line (37) 

X'(r)P(X) = r 2 — singular point 

In Section 4 we found the asymptotic expression for probability density function F Q = 
J F(z, fi,x)dfi ( |32l) . ( |33l) . Integrating it with respect to z, we obtain the function P(x): 

P(x) = J2p a x 3 - a a = A + 3n n = 1,2,3,... (38) 

a 

Combining (1381) with (1371) and integrating (j3"Tj) in the vicinity of the singular point, we get 
in the leading asymptotic term 

X(z) ~ z~t — singularity of vortex surface (39) 

_ 2 

X(r±) ~ r_!_ 3 — singularity of vortex line (40) 

X(r) ~ r _1 — singular point 

We see that in the case of singular point, X(r) diverges less than any singularity specified 
by the Laplace equation. This means that the singularity of such kind cannot exist as the 
isolated point. 

The most divergence in ([40]) is given by the singularity along the filament. We shall 

_ 2 

show that the dependence X(r±) ~ r ± 3 corresponds to the Kolmogorov law. Consider the 
correlator of transverse velocities 

K(r ± ) = <(v ± (r,t)-v ± (0,t)) 2 > (41) 
13 



Here is the velocity component perpendicular to the current line According to the 
definition (jlj), the expression Vj_(r, t) — vj_(0, t) is identically equal to w x r in the vicinity 
of the vorticity's singularity. Hence, as r — > we get 

K(r±) oc u 2 r\ oc r 2 / 3 (42) 

Since the direction of the line is arbitrary, the turbulent pulsations change it. From the 
isotropy it naturally follows that the space average over the main scale of the turbulence is 
equivalent to the average over the angles. The expression ( 14*21) then transforms to 

K(r) oc r 2/3 (43) 

So, the main expression (CQ), ( 1431) follows naturally from our consideration. It means that 
Kolmogorov's correlation function (14*31) is determined by the system of vortex filaments. The 
input of regular part of the velocity to the correlation function (oc r 2 ) is negligible. 



6 Conclusion 

In this paper we investigated the small-scale structure generation in the developed turbu- 
lence. Here we summarize and discuss briefly the main results. 

1. Vortex structures 

We considered the Navier-Stokes equation as v — > 0. We derived the equations describing 
the growth of small-scale pulsations along the Lagrangian trajectory under the action of the 
large-scale turbulence. We showed that the small-scale part of vorticity grew exponentially in 
time. This growth led to formation of a system of filaments and surfaces where the vorticity 
grew intensively. We derived the characteristic parameters of the vortex structure growth in 
time. 

2. Singularity 

We showed that in the non-dissipative limit v — > the absolute value of the vorticity uo 
tended to infinity along the vortex filaments as t — > oc. Notice that constructing the proba- 
bility density function (1201) we linearized the hydrodynamic equations near the Lagrangian 
trajectory ([9]). We obtained the exponential growth of vorticity which may be cut either by 
nonlinear corrections or by viscosity of the flow. Let us discuss now the nonlinear correc- 
tions. The feedback effect of the small-scale pulsations could be estimated by comparison of 
the energy density of small-scale pulsations with the energy of large-scale ones. The energy 
density of the main pulsations on the scale L is 

£o = \ P U 2 

here U is the velocity of the large-scale pulsations. From (1391) we estimate the velocity of 
the filament having width tq\ v n ~ U {j-q/L) 1 ^. Taking into account the part of volume 
occupied by the filaments, we obtain the relation between the energy densities: 
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Here Nj is the ratio of the number of filaments to the number of large-scale vortex in the 
volume unit. 

Thus, we see that the feedback effect of small-scale pulsations on the large-scale ones 
is insignificant. Hence, the singularity should be cut off by viscosity. The situation is 
quite analogous to that in supersonic hydrodynamics: the singularities (strong and weak 
discontinuities) in Euler flow are cut off by viscosity. 

3. Correlation function 

We found the solution of the equation describing the vorticity distribution in the vicinity 
of the vortex filaments. The solution had the form oc r 1 / 3 in the plane perpendicular to 
the vortex filament. This was the solution that determined the form of the pair correlation 
function in small scales. Thus, we found the velocity correlation function (CD), ( 1431) in the 
steady-state turbulent flow directly from the Navier-Stokes equation in non-viscous limit 
v -> 0. 

4. Intermittency 

According to the Kolmogorov-Richardson assumption the energy flux in turbulent flow 
cascades from larger scales to smaller ones and dissipates at the smallest scales uniformly in 
space and time. Landau pointed out that this assumption was controversial (see [3] §6.4). 
Gurvich in 1960 [18] . and later the other researchers discovered experimentally a very strong 
time and space inhomogeneity of velocity and energy flux. This property of turbulence is 
called intermittency. Variety of approaches to this effect was considered by many authors 
(see monographs [2], [3] and citations therein). 

Let us list the intermittency features that follow from the presented theory. 

1) The vorticity distribution in space is very inhomogeneous. Near the vortex axis it 
could possess the value many times exceeding its average. 

2) Even moments of the correlation functions should grow with number of the moment. 

3) The energy dissipation in the developed turbulent flow is localized near the axes of 
vortex filaments and the vortex surfaces. It is distributed very inhomogeneously in space 
and time due to the vortex structures motion. Besides, the strong nonuniform dissipation is 
the most pronounced manifestation of the intermittency [19J. 

Note that the filaments give the maximum degree of singularity as t — > oo and are respon- 
sible for the form of the pair correlation function. However, the surface-type singularities 
may affect the dissipation process, since they could occupy a significant part of the volume 
of the flow. 

The authors are grateful to V.L. Ginzburg for the attention to this work, and to A.S. 
Gurvich, V.S. Lvov, E.A. Kuznetsov, S.M. Apenko, V.V. Losyakov and M.O. Ptitsyn for 
useful discussions. 

This research was partially supported by the RAS Presidium Program "Mathematical 
methods in nonlinear dynamics" . 

Appendix 

Let us consider an axially symmetric flow. The hydrodynamic equations in the cylindric 



15 



coordinates take the form 



dv r dv r 
dt r dr 


dv r 
dz 


<p_ _ 

r 


dv 
dr 


(Al) 


dv^ | ^ dv^ 
dt ' " r dr 


, dv^ 
dz 


r 


= 


(A2) 


dv z dv z 
dt r dr 


dv z 
dz 


dp 

dz 




(A3) 


id. . 

r dr 


dz 






(A4) 



Here v r , v^, v z are the radial, asimutal and parallel to the cylinder's axis velocity components, 
respectively 

We search a solution of the system (A.l) - (A. 4) in the form 

v<f, = ur , v r = ar , v z — bz (^4-5) 
Then the pressure should be 

. PAt) 9 P 2 (t) 9 

p(r,z,t) = ^r 2 + ^z 2 

From (A. 4) follows a relation between a and b: 

2a + 6 = 0, (A6) 

This relation expresses the volume conservation in the liquid. Indeed, let us consider a 
cylindric drop with radius R(t) and length Z(t). Then from (A. 5) follows 

R = a(t)R, Z = b(t)Z, 

hence R(t)=R exp^J a(ti)dt^j , Z(t) = Z exp^J b{ti)dt^ 
The cylinder volume at arbitrary time t is 

nR{t) 2 Z{t) = 7r-RgZ„exp / (2a(*i) + 6(*i)) dt x = 7iR 2 Z . 

Jo 

We see that volume conserves. For example, if b > then the cylinder stretches, and its 
transversal radius decreases. 

Combining (A. 5) with (A.l) - (A. 4), we obtain a system of ordinary differential equations: 

a + a 2 — to 2 = —Pi 

u + 2aiu = (A7) 
b + b 2 = -P 2 
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Note that the system allows one arbitrary function of time. Actually, with account of 
(A. 6) we have four equations and five unknown functions: a , b , u , Pi , P-i- Without loss of 
generality one can choose Pi($) as such arbitrary function. We also note that the change 
of vorticity u{t) is connected unambiguously with change of the "cylinder length" Z(t): 

u(t) = U Z(t)/Z . 

Differentiating the second equation of the system (A. 7) and substituting other equations, 
we get 

u> = -P 2 {t)uj 

This equation is a particular case of (Q. We assume that P^if) is rather complicated "ran- 
dom" function and its time average is zero. Then the time intervals when P2{t) > and 
Pzit) < are equally probable. However, at P2(t) > the function uj(t) oscillates, the 
oscillation amplitude changing weakly. To the contrary, at P2{t) < the function u(t) 
grows exponentially. It is clear that in the average the value u grows. Since uo and Z are 
proportional, such growth means "stretching out" the cylinder. 
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